Skip to main content

Point Mass

Point Mass

Usage

point  
-0.03 #x-coord (= x centroid)
-0.04 #y-coord (= y centroid)
0.10 #theta_e2 (= $\theta_E^2$, strength of the point mass, log option available)

The deflection angle for a point mass at the origin is defined as:

α1=θE2r2x1andα2=θE2r2x2 \alpha_1 = \frac{\theta_E^2}{r^2} x_1 \quad \text{and} \quad \alpha_2 = \frac{\theta_E^2}{r^2} x_2

where:

r2=x12+x22r^2 = x_1^2 + x_2^2

Here, θE2\theta_E^2 is the strength of the point mass, related to the physical mass MM via:

θE2=4GMc2Dd\theta_E^2 = \frac{4GM}{c^2 D_d}

The factor (DdsDs)\left( \frac{D_{ds}}{D_s} \right) is incorporated separately in the config to convert this to a scaled deflection angle.

Calculating Mass MM

To compute the mass MM from a given θE2\theta_E^2, use:

θE2DdsDs=4GMc2(DdsDsDd)\theta_E^2 \cdot \frac{D_{ds}}{D_s} = \frac{4GM}{c^2} \cdot \left( \frac{D_{ds}}{D_s D_d} \right)

Thus, with known distances from cosmology, one can back out the physical point mass MM.

Einstein Radius

  • θE\theta_E is the Einstein radius of the point mass assuming a source at z=z = \infty.
  • To get the Einstein radius for a source at finite redshift:

θE,eff2=θE2DdsDs\theta_{E,\text{eff}}^2 = \theta_E^2 \cdot \frac{D_{ds}}{D_s}

Notes

  • The third parameter theta_e2 (i.e., θE2\theta_E^2) supports the log: option.
  • When using log: (log base 10):

θE2=10theta_e2_gleeconfig_value\theta_E^2 = 10^{\text{theta\_e2\_gleeconfig\_value}}

  • Units: θE\theta_E (and thus θE2\theta_E^2) should be in arcsec, to match the unit convention of other lens models in GLEE.
  • log: option is avialable for #theta_e2